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Abstract. We propose a mechanism for efficient heating of the solar chromosphere, 
based on non-ideal plasma effects. Three ingredients are needed for the work of this 
mechanism: (1) presence of neutral atoms; (2) presence of a non-potential magnetic 
field; (3) decrease of the collisional coupling of the plasma. Due to decrease of col- 
lisional coupling, a net relative motion appears between the neutral and ionized com- 
ponents, usually referred to as "ambipolar diffusion". This results in a significant en- 
hancement of current dissipation as compared to the classical MHD case. We propose 
that the current dissipation in this situation is able to provide enough energy to heat the 
chromosphere by several kK on the time scale of minutes, or even seconds. In this pa- 
per, we show that this energy supply might be sufficient to balance the radiative energy 
losses of the chromosphere. 



1. Introduction 



The degree of plasma ionization in the lower solar atmosphere - photosphe re and chro- 
mosp here - is very small. Using VAL-C model atmosphere as a reference dVernazza et all 
198 lh . the abundance of ionized atoms, relative to neutral atoms, is as low as 10~ 4 at 
heights of temperature minimum, and it remains always well below unity even at larger 
heights. In addition to that, the collisional coupling of the plasma becomes less im- 
portant with height. A simple calculation reveals that cyclotron frequency of hydrogen 
ions may exceed the collisional frequency already in the lower photosphere, for values 
of magnetic field strength expected for the quiet photosphere, see Figured] 

These two factors may lead to a break of the assumption underlying magnetohy- 
drodynamics (MHD) and lead to new effects, not taken into account in the classical 
approach, as ambipolar diffusion. In astrophysics, ambipolar (or neutral) diffusion usu- 
ally refers to the decoupling of neutral and charged components. Ambipolar diffusion 
causes the magnetic field to diffuse through neutral gas due to collisions between neu- 
trals and charged particles, the latter being frozen-in into the magnetic field. 

There is an increasing number of evidences for the importance of deviations from 

MHD in different situations. The presence of neutral atoms in partially ionized plasmas 

significantly affects wave excita tion and propagation ( Kumar & Roberts 20 031; iKhodachenko et al 
2004 120061 ; iForteza et al.l2007l : IVranies et alJl2008l ; ISoler et al.l2009Ll201d:IZaaarashvili et all 
2011 ) . It is also important for magnetic reconnection (Brandenbur g & Zweibellll994l 
"1 ISakai et al.ll2006l ; ISmith & Sakaill2008l ; ISakai & Smithl l2009h Non-ideal plasma 
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Figure 1. Estimation of the collisional and cyclotron frequencies for electrons (black 
lines) and hydrogen ions (red lines) in the photosphere and chromosphere of the Sun. The 
magnetic field strength is assumed to vary with height as B(z) = Boexp(-z/M B ), with 
B = 100 G and H B = 600 km. 



effect s can modify the equilibrium balance of photospheric flux tubes jK hodache nko & Zaitsev 



2002) and chromospheric structures, such as prominences (lArber et al.12 009: Gil bert et al 



2002). Another phenomenon potentially affe cted by non-idea l plasma effects is mag- 



netic flux emergence dLeake & Arbej|2006l : lArber et al.ll2007l) . Despite this increasing 
evidence, we are far from a complete understanding of th e influence of these effects. 

In this paper, we continue the investigation started in lKhomenko & Collados! (120121) . 
There, we studied the consequences of the ambipolar diffusion into the heating of the 
magnetized solar chromosphere. In the presence of neutrals, the ambipolar (or neutral) 
diffusion is orders of magnitude larger than the classical Ohmic diffusion, leading to 
efficient Joule dissipation of electric currents. Our calculations have demonstrated that 
just by existing relatively weak (10^10 G), non-force-free magnetic fields, the chromo- 
spheric layers above 1000 km can be efficiently heated by current dissipation reaching 
an increase of temperature of 1-2 kK in a time interval of minutes. The work of ambipo- 
lar diffusion would stop when all the atoms become ionized or when the magnetic field 
becomes force-free (/ || B). We proposed that this heating mec hanism may be efficient 
enoug h to balance the radiative losses of the chromosphere (see Khomenko & Collados 



20121 ). Here we explore this possibility by including the radiative damping in our initial 



calculations. 



2. Equations and numerical solution 



We solve numerically the quasi-MHD equati ons of conservation of mass , momentum, 
internal energy, and the induction equation (IKhomenko & Colladosll2012r) . 



| + ^(piO = 



(1) 



Du 
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where the following definitions are used: 



— = V x [(u x 3) - rjp f- TjAfJ-oJl] 



P= ^ Pa', it=- ^ (p ff 4); f=en e (jti-ti e ); p = ^ p ff (2) 



=n.i.e 



and 7 ± is the component of the current perpendicular to the magnetic field. These equa- 
tions are produced by summing up the equations for three different species (electrons 
(e), hydrogen ions (/) and neutral hydrogen («)). When deriving them we neglected the 
non-diagonal components of the pressure tensor and assumed that the diffusion veloci- 
ties w a = u - u a (a = e,i,n) are small, neglecting terms containing w 2 . In the Ohm's 
law we neglected the time variation of relative ion-neutral velocity the effects 

on the currents by partial pressure gradients of the three species, and the gravity force 
acting on electrons. The Hall term of the Ohm's law does not appear in the energy 
equation and, consequently, has no impact on the thermal evolution of the system. For 
consistency, we removed the Hall term from the induction equation as well. The Ohmic 
and ambipolar diffusion coefficients are equal to: 

m e (v ei + v en ) (Pn/p) 2 \B\ 2 

n = — 2 ; = t 1 r~ (3) 

e z n e p (piV in + PeV en )P0 

and the collisional frequencies are: 



3/2 
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nm in V nm m imje^ \ 2nk B T ) 

where m ; >, - m ; m„/(m/ + m n ) and m en - m e m n /(m e + m n ). The respective cross sections 
are cr in = 5 x 10~ 19 m 2 and <j m = 10~ 19 m 2 . A is the Coulomb logarithm. 

After subtractin g the equilibrium cond i tions, these equation are solved by means 
of our code mancha ( Khome nko et al.ll2008l : iFelipe et al.lEoTol : iKhomenko & Colladosl 



20121) with the inclusion of the physical ohmic and ambipolar diffusion terms in the 
equation of energy conservation and in the induction equation. As our code propa- 
gates non-linear perturbations to the magneto-static equilibrium, we treat the diffusion 
terms as perturbations. The radiative losses term <2rad in the energy equation calculated 
from the solution of Radiative Transfer Equation in Local Thermodynamic Equilibrium 
(LTE) and grey approximation for the opacity dependence on wavelength A. 

Grad = - f tfpA)dA; f A =[ htfWQ (5) 

J A JAn 

where Fa is radiative energy flux, l\ is specific intensity, p marks the direction and dQ. 
is a solid angle. Following the philosophy of our code, we also perturb the radiative 
cooling term, AQ m( \ = Q m( ± - <2 ra d(0), the zero-order term <2 ra d(0) being eliminated 
from the equations as a part of equilibrium condition. 
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Figure 2. Ionization fraction p e /p (left) and diffusion coefficients rj, tja in m 2 s -1 (right) as 
a function of height in the two flux tube models. Black thin lines: VAL-based tube; red thick 
lines: HOLMUL-based tube. 



The equations are solved in two spatial dimensions, though the vector quantities 
are allowed to have three dimensions (2.5D approximation). As the temperature and the 
ionization state vary with time, we recalculate the ionization balance of the atmosphere 
at each time step, assuming LTE (Saha equations). We then update the neutral fraction, 
p n /p, needed for the calculation of the ambipolar diffusion coefficient (Eq.|3]>. 



3. Flux tube model 

In quiet regions of the solar chromosphere, diagnostic tools base d on the Hanle ef- 



fect point to magnetic fie l d strengths of the order of ten s Gauss (ITrujillo Bueno et al. 
l2005l : ICenteno et al.l2010l : lstepan & Truji llo Bue noll2010l) . To simulate such non-active 



chromospheric conditions, we used a 2nd-order thin magnet ic flux tube model as ini- 
tial atmosphere ( Pneuman et alj[l986l : iKhomenko et al.ll2008l) . The model represents a 



horizontally infinite series of flux tubes that merge at some height in the chromosphere, 
preventing them from excessive opening with height. This magnetic field configuration 
is non-force-free. Here we used two flux tube models, the main difference between 
them is their (horizontally homogen eous) temperature str ucture. One has the verti- 



cal tempera ture structure of VAL-C (IVernazza et all 1 1981). and the other has that of 



HOLMUL (IHolweger & Muellerlll974l) . The magnetic field strength is similar in both 
cases, decreasing with height, from about 800 G in the photosphere to 35 G in the 
chromosphere. In the rest of the text, we will refer to these models as VAL-based and 
HOLMUL-based flux tubes. 

While the VAL-C model atmosphere includes the chromospheric temperature in- 
crease, the temperature in the HOLMUL model monotonically decreases with height 
dictated by the conditions of radiative equilibrium. This peculiarity determines the ion- 
ization fraction and value of the ambipolar diffusion coefficients in the models (see 
Fig. |2]). In the HOLMUL-based tube, the ionization fraction value does not exceed 
10~ 3 , even in the chromosphere. The values of the Ohmic diffusion coefficient 77 (see 
Eq. O are similar in both models. The ambipolar diffusion ija is orders of magnitude 
larger than rj already above 1 Mm height. Due to the much lower ionization fraction 
in the HOLMUL-based model, the value of tja there significantly exceeds that of the 
VAL-based model. Such values of rjA imply important current dissipation on very short 
time scales and a much quicker dissipation is expected in the cooler HOLMUL-based 
model. 
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Figure 3. Simulation snapshots after 800 s of evolution. Top left: VAL-C based flux tube 
with <2rad = 0; top right: same but Q ri & ± 0; bottom left: HOLMUL-based flux tube with 
Grad = 0; bottom right: same, but Q ml i t 0. The background color is temperature, the scale is 
the same for all panels. Vertical white lines are magnetic field lines. Arrows show the velocity 
field. For better visual comparison we keep unchanged the temperature structure from to 
0.42 Mm, though the variations are present in the simulations. 



4. Heating of small-scale flux tubes 

Here we describe four simulation runs: (i) VAL-based flux tube with <2 ra d = 0; (ii) 
VAL-based flux tube with g rad + 0; (iii) HOLMUL-based flux tube with Q mA = 0; 
(iv) HOLMUL-based flux tube with <2rad ^ 0. The flux tube models are initially in 
magneto-static equilibrium, obtained without considering the diffusion terms. Without 
external perturbation, they do not evolve. After introducing the perturbation in the form 
of diffusion terms, we perturb the initial magnetic field structure via the induction equa- 
tion. Then, as time evolves, this perturbation translates to the rest of the variables of the 
system and they start to change. Thanks to the Joule heating term in the energy equa- 
tion (r]AftoJ±, Eq. [T]) the magnetic energy is efficiently converted into thermal energy, 
producing heat. This heat is balanced by the radiative cooling term <2rad- 

Figure [3] shows snapshot from the four simulation runs 800 s after the introduction 
of the perturbation. At this time moment, the temperature has increased at the upper 
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Figure 4. Temperature as a function of height at horizontal position 0.6 Mm, see Fig. 3. 
Top left: VAL-C based flux tube with (2 rat | = 0; top right: same but <2, ac i # 0; bottom left: 
HOLMUL-based flux tube with gad = 0; bottom right: same, but Qm& ^ 0. Different lines 
are separated 1 sec in time at the upper two panels, and 0.5 sec in time at the lower two panels. 
Progressively more red colors indicate larger times till 800 sec since the start of the simulation. 



layers of the flux tubes in all simulations, but by a different amount. A more detailed 
view of the temperature behavior is provided in Figures @] and [5j Fig. @] gives the 
height dependence of the temperature at a fixed horizontal position inside the flux tube 
(X=0.6 Mm, close to flux tubes walls), for different time moments. Fig. [5] provides 
the time evolution of the temperature at a fixed point in the chromosphere, for all four 
simulations. 

In all cases, the most important heating is achieved at the upper part of the domain, 
close to the tube borders. This behavior is expected because the term responsible for 
the heating (tjaIUqJ 2 ) is orde rs of magnitude larger at these locations (see Fig. 9 in 



Khomenko & Collados 2012). 



In the simulations with = (left panels of Figs.OH]) the relative temperature 
increase, achieved after 800 sec, is significantly larger in the cooler HOLMUL model. 
In the VAL-based model, the temperature at 1.8 Mm reaches 8200 K after 800 sec of the 
simulation, which is ~1600 K above its initial value. In the HOLMUL-based model, it 
reaches 12000 K, i.e. -8200 K above its initial value. This disparity in the amount of 
heating is readily understood given the order-of-magnitude different rjA values (Fig.|2]). 
The temperature is visibly enhanced above 1 — 1.2 Mm height. At these heights the 
action of the ambipolar diffusion becomes important in our model. 

There is a large difference in the time scale of heating between the VAL-base 
and HOLMUL-based models. We can roughly define a time scale as L 2 /tja, setting 
a characteristic spatial scale of the system to L = 10 5 m. Figure [6] shows temporal 
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Figure 5. Time variations of temperature at X = 0.6, Z = 1.8 Mm in VAL-based flux tube 
(left) and HOLMUL-based flux tube (right). Dashed lines: Q mi = 0; solid lines: Q mi + 0. 
First few time steps are shown by bullets, separated in time by 1 sec (left) and 0.5 sec (right). 




Figure 6. Time variation of the ionization fraction p e /p (left); ambipolar diffusion coef- 
ficient r] A (middle); and characteristic time scale L 2 /rj A (right) in the VAL-based flux tube 
(black) and HOLMUL-based flux tube (red), in the simulations with <2 ra( i = 0. The values at 
taken at X = 0.6, Z = 1.8 Mm. Bullets are separated 0.5 sec in time. 



variations of the time scale in the VAL-based and HOLMUL-based models. It also 
shows variations of the ionization fraction and t]a- In the HOLMUL-based model the 
characteristic time scale is initially as low as 10 sec, producing an almost immediate 
increase of the temperature from 3700 K to 5700 K (Fig. \5\ bullets on the right panel). 
This temperature rise causes an increase of the ionization fraction from 10~ 3 to ~ 0.5, 
and a drop of the tja some three orders of magnitude at the first few seconds. After 
the initial rapid variation, the evolution becomes smoother, with characteristic times 
around 100 sec. Note that VAL-based model does not show such quick changes at the 
beginning of the simulation, since tja is initially much smaller than in the HOLMUL- 
based model. 

While the temperature constantly increases in the simulations with Q m & = 0, there 
is an oscillatory-like balance established in the simulations with Q m d & (Fig. 01 right 
panels; Fig. [5] solid lines). After an initial increase, we observe damped oscillations of 
temperature converging to some constant value. For example, at the location X = 0.6, 
Z = 1.8 Mm, the value of temperature to which the simulations converge is 7100 K 
(VAL-based, 500 K above the initial value), and 6000 K (HOLMUL-based, 2300 K 
above the initial value). Thus, keeping in mind all the approximations and simplifica- 
tions of our modeling, we conclude that the Joule heating and radiative cooling terms 
in the energy equation can balance each other. 
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5. Conclusions 

We have performed numerical simulations showing that the solar chromosphere can 
be effectively heated due to the Joule dissipation of electric currents, enhanced in the 
presence of neutral atoms (ambipolar diffusion). Our main conclusions are: 

• The amount of heating and its time scale depend on the initial temperature of 
chromospheric magnetic structures. In cooler regions the heating can act ex- 
tremely rapidly, reaching a temperature increase of 2 kK in few seconds time, 
while in the hotter regions the heating time scale is of the order of minutes. 

• The Joule heating by ambipolar diffusion may be able to balance radiative losses 
of the chromosphere. Our simulations show that after a period of damped oscil- 
lations, the temperature stabilizes at some constant value. In the particular case 
of the simulations considered here this value is 6 - 7 kK at 1.8 Mm. 
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